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LANDSAT 2 PROJECT MO. 289GD 
4th quarterly progress report 
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Title Water. Utilization, Evapotranspiration and Soil 
Moisture Monitoring in tlie SDuth East Region of 
South Australia. 

Assigned Investigation No . i 2896D 
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Author' s Name ; K.R. McCloy , t 

K, John Shepherd 
J.C. Killick, 

Reporting Date ; Dec, 2, 1976. (FTD + 14 months) 

TECHNIQUES 

1. Material Received 
70 mm Format NIL 

9 track, 1600 B.P.I. CCT. Frames 2342-23411 

and 2359-23351 received by 
the last quarterly report are 
being used. Frame 2198-23424 
of 8 Aug. 1975 was damaged in 
transit to Australia and has 
subsequently been returned to 
NASA. We have not yet received 
a replacement. 

2. Status of Project 


(i) Linear Classifier. The Linear Classifier 

developed by McCloy and mentioned in previous 
reports has now been used to classify Bool 
Lagoon at tliree dates, 29th November, 1972, 
30th December, 1975, and 16th January, 1976, 
and compared with aerial photography. These 
results with the dune survey have been offered 
for presentation at the Remote Sensing of 
Environment symposium In April 1977. A 
Comprehensive Summary of the proposed paper 
is included. 
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Detailed Landcover data has been collated into 
maps for the two dates, 30th December, 1975, and 
16th January, 1976. These maps are being 
digitised to correspond to the Landsat imagery 
at these dates. The data will be used to 
ascertain tlxe linearity of naturally occurring 
pasture types, for classification, and then 
testing the reliability of that classification. 
Initial work by McCloy suggests that the change 
in response across pasture is approximately 
linear, and related primarily to plant stage 
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of grov/th. 

(iii) Because of the late receipt of digital data 
it is requested that the project be allowed 
to continue for another two reporting periods 
qnd that a final report be produced at that 
time. It is anticipated that the next three 
months will be involved in finishing digitisa- 
tion of Landcovar, and extracting initial results 
from tliis. The second reporting period vrauld 
be used to classify test areas and analyse the 
reliability of the results. 

Ill ACCOMPLISHMENTS 


Successful operation of the improved Classifier 
routines. 

IV SIGNIFICANT RESULTS 


Nil this Quarter. 
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INTRODUCT ipN 

Most {|lp,It.’|]. clnsfilflor nlf’nrJ Llinifi nHsiimt* point rljninos, witli n dorint'd 
oiivolopc. Tills piipor duHiM'llnns Lho ilovo 1 opiiiont and d Ifu inifioa Htunjofial nJ, 
applications of a linear clnssifier. 

Conceptually, Ihc clnssifh't has a nnnilM'r of operator deslp.nated nodes or 
lliiiJtIng surl'at’O typos spoclflod liy iliolr mean covariance arrays. Thosu n nodes 
define a linear surface in m dimensional space. 

This linear surface has both dimension and orientation, lienee it can be 
considered to be a linear vector surface. Any pixel containing proportions of 
the various surfaces will lie on this vector, with proportions being defined by 
the position of the pixel with respect to the nodes. 

Tor classifying scanner imagery, the pixel must be assumed to not lie on 
the surface. The algorithm uses a constraint to f ind the footpoint on the surface 
that is "closest'' to the actual pixel. The closest Is in geometrical terms and 
under certain conditions the footpoint is near, but not the actual closest point. 
These conditions will only allow determination of the closest point in an 
involved manner, which is not considered to he justified. , 

The separation of tlie pixel from this footpoint can then be tested against 
statistical or geometric tliresliolds to determine whether the pixel con be classed 
as belonging, or not belonging, to the vector surface. 

Tl. procedure determine.^ the proportions of the various nodet; represented 
in the pixel, so tliat those pixels classified as belonging to the vector can he 
sub-classified according to their proportions of tlie various surface types. 

The derivation is based on independant types of surface cover similar to 
those used by Horwltz and lialepka in their article "Estimating the proportions 
of objects within a single resolution element of a niultlspcctral scanner"-, 
presented at the 7th Symposium in this series. To be maximally usefull, the 
classifier should be capable of classifying between the limits of an environ- 
mentally changing type of surface. The paper discusses a number of situations 
when tlie variations in radiance with variations In environmental condition, do 
display linear characteristics, and can he reliably classified with this technique. 

THEORETICAL-DERIVATION 


The power received by a scanner from an IFOV covered by a uniform surface 
of radiance is given by:- 

>■), ■ o ^ '^d + *■ 

If the IFOV contains n surfaces of spectral radiances 
the power received is given by:- 

l-’i = ■'"n -/"n^ t, da dO^' + N2i h, t da dO 

X o o 'A 0 d o o A 0 d 
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TliRORRTICAl^ nKRTVATIOW Gont,. 


Now = Ai = nj, etc.*. and the scanner response values ate J4ncarly 
U A 

related to i-ndiance and hence power received, so that; 


R^ = ai Rij^ -i* a?. Rz^ 4-. . . .-i- a^ Riij^ 


. . .3 


R^ is the pixel response, tlie response of surface 1, etc.. 

Xf tliere are m wavebands W,X,Y, ...M and n surfaces l,2,3,...n then we 
would have the in equations 

R =^aiRi +aaRz +....+ a Rn 
w *w . *w n w 


R = ai Ri +32 Rz + + a Rn 

X X X n X 


...A 


R - ai Ri + az Ra +....+ a Rn 
m m m n m 

Now for any pixel, not on the vector, we want the footpolnt on the vector. 
The footpolnt will be that point for which 

ai + az + a 3 +. . . .q- a = 1 .• ... 5 

Substitute this into equation A to give: 

R - Rn = Hi (Ri - Rn ) d- a? (Rz " Rn ) q-....+ a - Rn ) 

w w n w w w n-j w w 

R - Rn - 3] (Ri - Rn ) + az (Rz ~ Rn ) . . .+ a (Rn-i - Rn ) 

X X XX XX n“l X X 
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and as long as M ^ (n-I) then these equations can be solved for aj, az,...an-i 
and by equation 5, for ap- 

(R ,R ,R ,...R ) are the actual pixel response values, 
w’ X y’ in' 

Generally there will be less than (M+1) nodes so tliat the solution needs to 
be done by least squares. Equation 6 can be stated in matrix notation as: 


D = 0 * U wliere U is the column vector of unkno\m proportions. 


T T 

0 *0*U = 0 *I) 

T 0 

Let 0 * 13 = A 

U = a“^ * o"'''’ * D 
Let a"* * o’*’ = B 
U = n * D 
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THEORETICAT., DERIVATION Cont.. 


Army B is constant for any vector and can be determined prior to classification, 
so that classification only depends upon the matrix multiplication and (n-1) 
subtractions to determine 3^,. The footpoint response values ore then determined 
by substituting the proportion values into equation 4. 

Often the footpoint will fall outside of the domain of the nodes, in which 
case one or more proportions will be negative. One or more may be positive. 

Tlie author decided to restrict the footpoint to within the domain of the nodes. 

This is done by making all negative proportions zero and dividing all new 
proportions by the sum of the new proportions. This approach has advantages 
and disadvantages, which will be discussed, the author believing that the advantages 
outweigh the disadvantages. 

v; 

If the covariance arrays for the nodes are given as COVi, COVj, , . ,C0V;i) then 
it can be shown, by the law of propagation of error, thatt 

2.2 2 
GOV (I,J) = a. COVi (I,J) + 32 COV 2 (I,J) a„ COVn 

for all (I,J). This gives the covariance array of the selected footpoint, so 
that the likelihood of the actual pixel belonging to the vector can be tested 
statistically, in a variety of ways. 

The statistical envelope defined in this way are parabolic surfaces on 
either side of the vector surface. The author believes that in many cases it 
is not convenient to either determine or estimate the covariance arrays and 
believes that testing the geometric distance between pixel and footpoint will 
not introduce significant errors, for a saving in classification time. 

ENVIRONMENTAb SUPPORT POR CEASSIPIER 


The turbid lakes of Western Victoria display linear relationships between 
the response bands through the series of lakes sampled. The response values 
are related to turbidity levels, although not necessarily in a linear manner. 

Profiles across a large swamp lagoon, Bool Lagoon, give a linear relationship 
between the bands, from water surfaces through to the dry reed beds. 

. RESULTS ACHIEVED 


The classifier had been ucacl in two studies, Bool Lagoon and classifying 
the Sand Dunes from coastal scrub near Lake Bonney on the South Australian Coast. 
Both exercises were successful! and will be reported in more detail at the 
symposium. The Bool Lagoon work was done at three dates, 29th November, 1972, 
30th December, 1975 and 16th January, 1976. 


